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We study the properties of the g-state frustrated bond percolation model by a Monte Carlo "bond 
flip" dynamics, using an algorithm originally devised by Sweeny and suitably modified to treat 
the presence of frustration. For q = 2 the model gives the cluster formulation of the Edwards 
Anderson spin glass. We analyze the percolation transition of the model, and find that it falls in 
the universality class of the q/2-state ferromagnetic Potts model. We then investigate the bond flip 
dynamics of the model, and find that, while for temperatures higher than the percolation transition 
T p the relaxation functions are fitted by a single exponential, for T < T p they show a two step 
decay, reminiscent of the relaxation of glass forming liquids. The long time decay (a-relaxation) is 
well fitted for T < T p by a stretched exponential function, showing that in this model the relevant 
mechanism for the appearing of stretched exponentials is the percolation transition. At very low 
temperatures the relaxation functions develop a long plateau, as observed in glass forming liquids. 



I. INTRODUCTION 

Cluster concepts have been extremely useful in crit- 
ical phenomena to elucidate the mechanism underlying 
a thermodynamical transition, by providing a geometri- 
cal interpretation of thermodynamic correlations. In the 
Ising model defined by the Hamiltonian 
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a correct definition of clusters was given by Kasteleyn 
and Fortuin Q and by Coniglio and Klein Q|. In the 
Coniglio-Klein approach, one throws a bond between 
pairs of parallel nearest neighbor spins with a probabil- 
ity p = 1 — e~ 2l3J , where (3 — l/fc^T. By summing over 
the spin configurations with the Boltzmann weight, the 
partition function of the model defined by Eq. (Q) can 
be written as a sum over bond configurations C, 
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where q — 2, b(C) is the number of bonds in the config- 
uration C, N(C) the number of connected clusters, and 
fi = log(j-^-) = log(e 2 ^ J — 1) is the chemical potential 
of the bonds. Thermodynamic averages can be related 
in this approach to corresponding percolative quantities. 
One finds that the clusters represent spin fluctuations, 
and percolate at the Ising critical temperature with Ising 
critical exponents. 

This approach can be extended to the Potts model, in 
which spins can have q ^ 2 different states. In this case 
the parameter q in Eq. (^) can assume a value q ^ 2, and 
for every value of q the percolation model has the same 
critical temperature and exponents of the corresponding 
Potts model. For q = 1 one recovers the random bond 



percolation. Sweeny ||] studied the weighted percolation 
problem defined by Eq. (||) by Monte Carlo techniques on 
a two-dimensional square lattice. He showed that a simu- 
lation based on this approach does not suffer from critical 
slowing down for q < 4. In a few Monte Carlo steps one 
can equilibrate even a very large lattice at the critical 
temperature. Thus he extracted informations about the 
critical point of the Potts model by measuring geomet- 
ric quantities like the mean cluster size at the transition 
point. Later, the cluster approach was further elabo- 
rated by Swendsen and Wang Q and Wolff || , by imple- 
menting an efficient cluster dynamics to simulate directly 
the spin model (Q), which drastically reduces the critical 
slowing down of conventional spin flip simulations. 

The cluster approach has been extended also to frus- 
trated systems like spin glasses HJtJ . With respect to the 
ferromagnetic case, new features appear, due to the phe- 
nomenon of frustration. In fact the percolation model 
that one obtains has the same complexity of the original 
spin model, and is not useful to define fast Monte Carlo 
dynamics, as in the unfrustrated case. Nevertheless, it 
represents an interesting tool to investigate the proper- 
ties of frustrated spin systems from the geometrical point 
of view. Moreover, it can be considered on its own as a 
model of percolation in a frustrated medium, that makes 
it of relevance in the study of systems in which frustra- 
tion and connectivity play a central role, as the structural 
glasses. 

To illustrate the cluster approach to frustrated spin 
models, let us consider the Ising spin glass, defined by 
the Hamiltonian 
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where are quenched random variables that can have 
the values = ±1. As in the ferromagnetic case, one 
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throws a bond between nearest neighbor spins that sat- 
isfy the interaction (in this case if tijSiSj — 1) with 
probability p = 1 — e~ 2 @ J . The crucial difference with 
the ferromagnetic case is that in general not all the inter- 
actions can be satisfied simultaneously. Indeed a closed 
loop such that the product of along the loop is nega- 
tive does not admit any spin configuration satisfying all 
the interactions, and is called frustrated. Since bonds 
can be put only between spins that satisfy the inter- 
action, one cannot put bonds on the lattice that close 
a frustrated loop. In terms of bond configurations, the 
partition function of the Ising spin glass turns out to be 
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where q = 2, fj, = log(^-) = log(e 2 ' 3J - 1), b{C) is the 
number of bonds, and N(C) the number of clusters in 
the configuration C. Here the sum J^c * s extended to all 
the bond configurations that do not contain a frustrated 
loop. Therefore, in this cluster formalism, the only differ- 
ence between the Ising spin glass and the ferromagnetic 
model is that, in the ferromagnetic case (^|) the sum is 
over all the bond configurations, while in the spin glass 
(f|) due to the geometrical constraint the sum is restricted 
to a subset of the bond configurations. In particular the 
ground state at T = (// = oo) is obtained by maximiz- 
ing the number of bonds, under the constraint that the 
bond configuration does not contain a frustrated loop. 

Due to frustration, clusters defined in the spin glass 
model no longer correspond to thermodynamical fluctua- 
tions. In fact the correlations between spins can be either 
positive (if they propagate along a path that contains an 
even number of negative interactions) or negative (if the 
path contains an odd number of negative interactions), 
so they interphere and are canceled out at least in part 
On the other hand connectivity is always positive 
and, as a result, clusters in the spin glass model perco- 
late at a higher temperature T p respect to the critical 
temperature T$g- I n the three dimensional Ising spin 
glass simulations show that T p ~ 3.95J/ks IjJ, while 
T SG ^ l.HJ/fcfl §. 

As in the ferromagnetic case, the model defined by Eq. 
(|J) can be extended to values q ^ 2 of the spin multiplic- 
ity. We call this model the "q-state frustrated percolation 
model". For q integer and even, Eq. (^) is the partition 
function of the model defined by the Hamiltonian 10 
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where <7j = 1 ... q/ 2 are Potts spins, and S aiCrj is the Kro- 
neker's delta. In analogy to the case 9 = 2, the model 
is expected to have two transitions, one at a tempera- 
ture Tsg(q) corresponding to freezing of Ising spins, and 
the other at a temperature T p (q) > Tsg(q) correspond- 
ing to the percolation of the clusters. A renormalization 
group calculation carried over a hierarchical lattice flO] 



has confirmed this expectation, and has shown that the 
transition at T$g should be in the universality class of 
the Ising spin glass, no matter what is the value of q, 
while the percolation transition at T p should be in the 
universality class of the q/2 state Potts model. For q = 1 
the model assume a very simple form, as the factor q N ( c > 
disappears. The resulting model has been called "frus- 
trated percolation" , and can be viewed as a simple model 
of percolation in a frustrated medium. Despite its sim- 
plicity, its dynamical properties exhibit a complex behav- 
ior, with features in common with both structural glasses 
and spin glasses. 

In this paper, we perform a Monte Carlo study of the q- 
state frustrated percolation model on a two-dimensional 
square lattice, using a "bond flip" dynamics, in which 
bonds are added and removed from the lattice with ap- 
propriate probabilities, in order to satisfy the principle 
of detailed balance. To do this, we have realized an algo- 
rithm that allows to determine the connectedness of two 
given sites, and the presence of frustrated loops, with a 
time that in the worst case (at the percolation tempera- 
ture) scales only with the logarithm of the lattice size, as 
we describe in Sect. 0. The algorithm is a modification 
of the algorithm used by Sweeny in the ferromagnetic 
case M. In Sect. |jll] we report our results on the perco- 
lation transition of the model for q = 1, 2,4. Note that 
while for q multiple of 2 the percolation transition can 
be studied also by conventional spin flip [^2[ , for other 
values of q the "bond flip" dynamics is the only way to 
simulate the model. In Sect. |y| and fv| we study the dy- 
namical properties of the model, analyzing the relaxation 
functions of the number of bonds. 



II. MONTE CARLO ALGORITHM 

We have implemented a Monte Carlo algorithm to sim- 
ulate the bond percolation model defined by Eq. (Q) on a 
two-dimensional square lattice, which can be applied for 
any value of the parameter q € [0, oo). The interactions 
6ij between pairs (ij) of nearest neighbor sites are set at 
the beginning to a value +1 or —I randomly, with equal 
probability. These variables are quenched, and their state 
is not changed by the dynamics. 

Each edge of the lattice, that is each pair of nearest 
neighbor sites, can be in two possible states: connected 
by a bond or not. At each step of the dynamics, we 
try to flip the state of an edge chosen randomly, with 
a probability determined in such a way to satisfy the 
principle of detailed balance. If we try to remove a bond 
from the system, or to add a bond that do not close a 
frustrated loop, then the probability of carrying out the 
"bond flip" will be 
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where 8b is the change in the number of bonds, 5N is the 
change in the number of connected clusters. If we are 
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trying to add a bond that closes a frustrated loop, then 
we have simply P flip = 0. A Monte Carlo Step (MCS) is 
defined as 2V single bond flip trials, where V = L 2 is the 
total number of sites, and 2V the total number of edges 
of the lattice. 

The nontrivial point here is to determine the change 
in number of connected clusters, and to verify if a bond 
added between two given sites closes or not a frustrated 
loop. To do this, we have used the algorithm used by 
Sweeny in the ferromagnetic case ||, suitably modified 
to treat the frustration occurrence. Consider a two- 
dimensional square lattice, together with its dual lattice. 
If a bond is present on the lattice, then its dual bond 
is absent, and viceversa. The boundaries between con- 
nected clusters on the lattice and on its dual will form a 
collection of closed loops, as shown in Fig. 0(a). These 
loops are represented in the computer as chains of point- 
ers. Each site on the lattice has four pointers adjacent to 
it, as shown in Fig. |l](b). At the beginning of the sim- 
ulation, pointers are organized in a hierarchical way, by 
giving them a defined "level" , that do not change in the 
following. A fraction (4~™ — 4~(™+ 1 )) of the pointers are 
at level n = 0, 1, . . . , n max — 1, and 4^"™x a ^ i eve i n max , 
where level n max is chosen so that it counts not more than 
four pointers. Chains then are formed by making each 
pointer point to other two pointers of the chain, one in 
the direction of the arrows, called the ALONG pointer, 
and one in the opposite direction, called the UP pointer. 
The ALONG pointer must be at least at the same level 
of the one pointing to it, and the UP pointer at a higher 
level (except if there are no higher level pointers in the 
chain) , so they in general do not correspond to the near- 
est pointers in the chain. 

When we add a bond to the lattice, (and remove its 
dual), two things can happen: the bond links two sites al- 
ready belonging to the same connected cluster, and there- 
fore 5N = 0; the bond links two previously disconnected 
clusters, and S N = — 1. In the first case the bond will 
cut a single chain into two distinct chains, see Fig. ||(a), 
while in the second case it will join two distinct chains 
into a single chain, see Fig. |(b). When we remove a 
bond from the lattice, (and add its dual), it happens the 
other way round. Using this auxiliary data structure, one 
can determine if two given chain pointers belong to the 
same chain or not, and cut and rejoin chain segments, in 
a CPU time that grows only with the logarithm of the 
chain length. 

To simulate the frustrated percolation model, we must 
also determine if a bond added to the lattice closes a 
frustrated loop or not. To do this, we must be able to 
count the number of antifcrromagnetic bonds encoun- 
tered along a path that joins two given sites A and B. 
This can be done if every chain pointer contains infor- 
mation about the number of antiferromagnetic bonds 
"skirted" when one traverses the chain to its UP pointer. 
Call this number the "phase" of the pointer respect to its 
UP pointer. We then go on jumping from each pointer 
to its UP pointer, and adding the relative phases, until 



we reach a reference pointer R in the chain. Then the 
number we seek for is found as the difference between 
the phases of two pointers adjacent to the sites A and 
B, respect to the reference pointer R, as shown in Fig. 
||. This reference pointer is chosen between the highest 
level pointers in the chain, and is the UP pointer of all 
the pointers belonging to the highest level in the chain. 
When we cut and rejoin chain segments, we must up- 
date coherently the relative phases of the chain pointers 
involved, and assure that there is one and only one ref- 
erence pointer per chain. 

III. THE PERCOLATION TRANSITION 

We have studied the percolation transition of the 
model defined by Eq. ([|), for q = 1,2,4. For each 
value of q we have simulated the model for lattice sizes 
L = 32, 64, 128. The histogram method JTJj was used to 
analyze the data. In this Section, we use the probability 
p as the independent variable. It is connected to the tem- 
perature via the simple relation p = 1 — e~ 2 @ J . For each 
value of q and L, sixteen probabilities were simulated 
around the percolation transition point, taking 10 3 MCS 
for thermalization, and between 10 4 MCS (for L = 128) 
and 10 5 MCS (for L = 32) for the acquisition of the his- 
tograms. Histograms were taken of the number of bonds; 

of the mean cluster size, defined as — s 2 n Sl where n s 

s 

is the number of clusters of size s; and of the occurrence 
of a spanning cluster, defined as a cluster that spans from 
the bottom to the top of the lattice. These histograms 
were used to calculate the average number of bonds (6) , 
the fluctuation in the number of bonds (b 2 ) — (b) 2 , the 
average mean cluster size x> and the spanning probabil- 
ity -Poo, in a whole interval of probabilities around the 
percolation transition probability. 

The spanning probability Pca(p) and the mean clus- 
ter size x{p)i as a function of the lattice size L and of 
the probability p = 1 — e~ 2 @ J , around the percolation 
probability p c should obey the scaling laws [l^Jl5|] 

Poo(p) siP o[L 1/v (p-Pc)], (7a) 

X (p)^L^x[£ l7 >-Pc)], (7b) 

where 7 and v are the critical exponents of mean clus- 
ter size and connectivity length, P^ and \ are universal 
functions. Given Eq. (|7a|), the value of the transition 
probability p c can be evaluated from the point at which 
the curves Poo(p) for different values of L intersect. In 
Fig. ||(a) are plotted the measured curves -Pqo(p) for 
q = 1 and L = 32, 64, 128. One must extrapolate the 
value at which curves for L, L' — > 00 intersect. Then 
we have evaluated the critical exponent 1/v, by choosing 
the value that gives the best data collapse of the curves, 
when one plots Pooijp) in function of L x ^{p — p c ), see 
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Fig. 0(b) . The errors on p c and 1/v were evaluated as 
the amplitudes of the intervals within which a good data 
collapse was obtained. Given Eq. (|Tb|), the mean clus- 
ter size x(p) nas a maximum that scales as L 1 ^, see Fig. 
~\(a,). From a log- log plot of ^ max in function of L, we can 
extract the exponent "f/v, making a linear fit, as shown 
in Fig. 0(b). Results are summarized in Tab. |. These 
results are in good agreement with the prediction that 
the percolation transition of the model falls in the uni- 
versality class of the g/2-state ferromagnetic Potts model 
jnj , whose critical exponents are reported in Tab. || . 

Renormalization group calculations show that the 
model should also have a singularity in the free energy 
density F(j>) at p c (q), with a singular part F s - mg (p) ~ 
A{q){p — p c ) 2 ~ a ( q / 2 \ where a(q/2) is the specific heat 
exponent of the g/2-statc ferromagnetic Potts model. In 
particular the model with q — 4 should have a singularity 
corresponding to the ferromagnetic Ising model, that is 
a logarithmic divergence of the second derivative of the 
free energy (specific heat). From Eq. (0) it is possible 
to show, that the derivatives of —(3F — V~ l \ogZ re- 
spect to fj, are equal to the cumulants of the distribution 
of the number of bonds, divided by the total number of 
sites V. As /i is a regular function of the probability p 
for < p < 1, we conclude that the model with q = 4 
should have a divergence in the second cumulant, that 
is the fluctuation, of the number of bonds divided by V. 
For finite size systems, we expect to see a peak whose 
maximum scales as log(L). In Fig. |^ we show our re- 
sults for q — 4 and L — 32, 64, 128. It is evident that 
there is a divergence, but statistical errors do not allow 
to distinguish between logarithmic and a weak power law 
divergence. 

For q = 2 one would expect a singularity in the free 
energy, at the percolation transition, characterized by an 
exponent 2 — a(l) = 8/3, that is a divergence in the third 
cumulant of the number of bonds. For the Ising spin 
glass (q — 2) this would imply a divergence in the third 
cumulant of the distribution of the energy. However sev- 
eral arguments, including a renormalization group calcu- 
lation |LlJ , predict that the singularity could be canceled 
out by the vanishing of the prefactor A(q) for q — > 2. 
This prediction could be checked by showing that the 
third cumulant does not diverge, but much more exten- 
sive simulations are needed to verify this prediction. 



IV. ONSET OF STRETCHED EXPONENTIALS 

We have studied the dynamical properties of the q- 
state frustrated percolation model for q = 1, 2, and 4, by 
calculating the autocorrelation function of the number of 
bonds. This is defined as 
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dimensional square lattice of size 32 x 32, and took 10 5 
MCS for thermalization, and between 10 5 and 2 x 10 6 
MCS for acquisition. All the functions were then aver- 
aged over 16 different configurations of the interactions 
Cij. Errors are evaluated as mean standard deviation of 
this last averaging. 

In Fig. 0(a), ||(a), and 0(a), we show respectively the 
results for q = 1, 2, and 4. For all the three values of the 
multiplicity q considered, we observe the same behavior. 
For high temperatures, the functions show a single expo- 
nential decay, while at lower temperatures they show in- 
stead a two step decay, reminiscent of what one observes 
in glass forming liquids at low temperature. This behav- 
ior can be explained by the presence, at temperatures 
below the percolation transition, of a rough landscape of 
the free energy in configuration space, with many minima 
separated by high barriers. The short time decay corre- 
sponds to the relaxation inside the single valley, while the 
long time tail is due to the tunneling through barriers and 
final decay to equilibrium (a relaxation). Furthermore, 
the second step of relaxation functions is well fitted by a 
stretched exponential function 



F(t) oc exp — 
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and is normalized so that F(0) = 1, while for t — > oo 
it relaxes to zero. We simulated the model on a two- 



where [3 is an exponent lower than one. 

In disordered and frustrated spin systems like spin 
glasses, simulated by conventional spin flip, the appear- 
ing of non exponential relaxation is believed to be caused 
by the existence of unfrustrated ferromagnetic clusters of 
interactions, see Randeria et al. ]l7] . Below the ferromag- 
netic transition temperature T c of the pure model, each 
unfrustrated cluster relaxes with a time that depends 
from its size. Due to the disorder of the interactions, 
the sizes of the unfrustrated clusters are distributed in 
a wide range, giving rise to a wide distribution of relax- 
ation times in the model. Therefore, according to this 
picture, the temperature T c of the ferromagnetic transi- 
tion of the pure model marks the onset of non exponential 
relaxation. 

The g-state frustrated percolation model is equivalent, 
for what concerns static properties, to a disordered spin 
system. In particular there are unfrustrated clusters of 
interactions with different sizes, due to the disorder of 
the variables ey. On the other hand, the bond flip dy- 
namics does not suffer from critical slowing down near 
the ferromagnetic critical point || , so we expect that the 
Randeria mechanism does not apply in this case. We 
have verified this point plotting the stretching exponent 
(3{T) as a function of the temperature T, for q = 1, 2, and 
4, as shown in Fig. 0(b), ||(b), and||(b), respectively. In 
the case of q = 2 and 4, it is quite evident from the data 
that the transition point of the pure model T c does not 
mark any change in the behavior of relaxation functions. 
Instead the temperature T p , at which clusters of bond 
percolate, appears as the point that marks the onset of 
stretched exponential relaxation. The case q — 1 was not 
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so evident from our results. For this reason we made a 
single simulation for a much larger system (L = 100), at a 
temperature T = 2.3 slightly higher than the percolation 
threshold (T p = 2.25). By fitting the function for times 
greater than 1.5 MCS, where F(t) < 0.05, we obtain a 
stretching exponent (3 — 0.97, definitely higher than that 
obtained at the same temperature for L = 32 (f3 = 0.86). 
Thus it seems that for q = 1 finite size effects at the 
percolation transition are more important, but neverthe- 
less the relaxation is asyntotically purely exponential for 
T > T p . In conclusion our results show that, for all the 
values of q studied, in this model the relevant mecha- 
nism for the appearing of non exponential relaxation is 
the percolation transition. 

V. DYNAMICAL PROPERTIES AT LOW 
TEMPERATURE 

We have evaluated the relaxation functions, in the 
model with q = 1 and 2, for very low temperatures. 
In Fig. |l(](a) and (b) we show the results for L = 40, 
q = 1 and 2, and for different temperatures, averaged 
over 16 different configurations of interactions. For the 
lowest temperatures, the functions do not relax smoothly 
to zero. This is clearly an effect due to the relaxation 
time being greater than the total time of the run, that 
was between 2 x 10 6 and 5 x 10 6 MCS. For such very low 
temperatures, the relaxation functions show a behavior 
very similar to what can be observed in glass forming liq- 
uids near or below the mode coupling theory transition 
temperature Jl8| |. A first short time decay is followed 
by a very long plateau, and eventually there is a final 
relaxation to equilibrium, for very long times. 

We compare this behavior to what one observes in the 
Ising spin glass simulated by conventional spin flip. In 
Fig. [ll] the relaxation functions of the energy, for the 
Ising spin glass on a two-dimensional square lattice with 
L = 40, and for the same temperatures of Fig. p"o|(b) , 
are shown. Note that in this case there is not a clear 
separation between the first short time decay and the 
long time tail of the functions, and the functions do not 
show any plateau. 



We have then studied the dynamical properties of the 
model. Differently from what happens in spin glass sys- 
tems, simulated by conventional spin flip, the transition 
temperature T c of the pure model does not play here any 
role in determining the dynamical behavior. Instead, the 
percolation temperature T p appears to mark the onset of 
two step decay, and stretched exponentials in autocorre- 
lation functions. 

At very low temperatures the autocorrelation functions 
develop a long plateau, as observed in glass forming liq- 
uids, and predicted by the mode coupling theory. This is 
a feature of the bond flip dynamics we have performed, 
while the spin flip dynamics of the spin glass model, that 
is thermodynamically equivalent to our model for q = 2, 
is very different, and does not show any plateau. 

This results show that the frustrated percolation model 
bridges the spin glass to the glass forming liquids. On 
one hand it is equivalent thermodynamically (for q = 2) 
to the Ising spin glass. On the other hand the model 
describes a bond packing problem, where the bond are 
subjected to the constraint that no frustrated loop can be 
closed. This makes the dynamics similar to that of glass 
forming liquids, where irregular molecules (or groups of 
molecules) move under the constraint of some kind of 
steric hindrance. 

To make more contact with glass forming liquids, the 
site version of the frustrated percolation model has been 
developed |]LS}] . In this model particles are allowed to 
diffuse on the lattice, under the constraint that no frus- 
trated loop can be fully occupied. Consequently it is 
possible to calculate the mean square displacement and 
the diffusion coefficient of the particles. Numerically it is 
found that also these quantities reproduce qualitatively 
the corresponding quantities measured in glass forming 
liquids. 
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VI. CONCLUSIONS 



We have studied the q-state frustrated percolation 
model, by means of a "bond flip" Monte Carlo dynamics. 
The model is equivalent from the thermodynamic point 
of view to the Potts spin glass model (|J) , which for q = 2 
coincides with the Ising spin glass. We have studied the 
percolation transition, that happens at a temperature T p 
greater than the spin glass temperature Tsg, an d found 
that it belongs to the universality class of the q/2-state 
ferromagnetic Potts model. 
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FIG. 2. Addition of a bond to the lattice: (a) the bond 
links two sites already belonging to the same cluster, and 
cuts a single chain into two distinct chains; (b) the bond links 
two disconnected clusters, and joins two distinct chains into 
a single chain. 
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FIG. 1. (a) Bonds, dual bonds, and chains. Open circles 
represent sites of the original lattice; solid squares sites of the 
dual lattice; solid lines bonds and dual bonds; dashed lines 
represent chains, (b) Pointers forming chains on the lattice. 
Each pointer points to an ALONG pointer in the direction 
of the arrow, and to an UP pointer in the opposite direction 
(not necessarily the nearest ones). 
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FIG. 3. Determination of the number of antiferromagnetic 
bonds along a path from site A to site B. Straight lines repre- 
sent ferromagnetic bonds, wavy lines antiferromagnetic ones. 
The number of antiferromagnetic bonds skirted traversing the 
chain in the UP direction (opposite to the arrows) from the 
pointer G adjacent to A to the reference pointer R is 3, while 
from H to R is 2. Making the difference, we find the number 
we seek, that is 1. 
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FIG. 5. (a) Mean cluster size x^(p) as a function of proba- 
bility p, for lattice sizes L = 32, 64, 128. (b) Maximum x™ ax 
of X^(p) as a function of lattice size L. 



FIG. 4. Measure of percolative quantities on the model 
with q = 1 and lattice sizes L = 32, 64, 128. (a) Spanning 
probability -Poo(p) as a function of probability p. (b) Poo(p) 
as a function of (p — p c )L 1//l/ , with p c = 0.589 and 1/v = 0.56. 




FIG. 6. Fluctuation (6 2 ) — (b) 2 of the number of bonds, 
divided by the number of sites V, as a function of p, for q — 4 
and lattice sizes (from bottom to top) L = 32, 64, 128. 
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FIG. 7. (a) Relaxation functions of the number of bonds 
in the model with L = 32, q = 1, for temperatures (from left 
to right) T = 2.7, 2.3, 2.0, 1.8, 1.5, 1.2. Solid lines are the 
stretched exponential fit functions, (b) Stretching exponent 
/3(T) as a function of temperature for L — 32 (stars). The 
open circle represents a single simulation made for L = 100 
and T = 2.3. Arrows mark the percolation transition T p and 
the critic temperature of the pure model T c . 



FIG. 8. (a) Relaxation functions of the number of bonds 
in the model with L — 32, q = 2, for temperatures (from left 
to right) T = 2.2, 2.0, 1.8, 1.5, 1.2, 1.0. Solid lines are the 
stretched exponential fit functions, (b) Stretching exponent 
(3{T) as a function of temperature. 
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FIG. 9. (a) Relaxation functions of the number of bonds 
in the model with L = 32, q — 4, for temperatures (from left 
to right) T = 1.9, 1.7, 1.5, 1.3, 1.2, 1.0. Solid lines are the 
stretched exponential fit functions, (b) Stretching exponent 
P(T) as a function of temperature. 



FIG. 10. (a) Relaxation functions of the number of bonds 
in the model with L = 40, q = 1, for temperatures (from 
bottom to top) T = 1.5, 1.0, 0.6, 0.5. (b) Model with L = 40, 
q = 2, for temperatures (from bottom to top) T = 0.6, 0.5, 
0.4, 0.35. 
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FIG. 11. Relaxation functions of the energy in the spin 
glass model with L — 40, simulated by spin flip, for temper- 
atures (from bottom to top) T = 0.6, 0.5, 0.4. 



q 


T P 


l/i/ 




l 


2.25 ±0.03 


0.56 ±0.01 


3.19 ±0.05 


2 


1.814 ±0.024 


0.77 ±0.02 


2.28 ± 0.06 


4 


1.47 ±0.04 


0.95 ± 0.04 


1.79 ±0.08 



TABLE I. Measured critical temperature and exponents 
for the percolation transition in the g-state frustrated perco- 
lation model. Temperatures are measured in unities of J/fcs. 



q 


T c 


1/u 


7 


OL 


1/2 


3.740 


0.5611 


3.2696 


— 1.5645 


1 


2.885 


3/4 


43/18 


— 2/3 


2 


2.269 


1 


7/4 


(log) 


4 


1.820 


3/2 


7/6 


2/3 



TABLE II. Critical temperature and exponents of the 
g-state ferromagnetic Potts model. Temperatures are mea- 
sured in unities of J /kn- 
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